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Abstract. We discuss various models and scenarios proposed to explain the nature of the V838 Mon type eruptions. 
In this class of eruptive objects we include: M31 RV (erupted in 1988), V4332 Sgr (erupted in 1994) and V838 
Mon (erupted in 2002). We concentrate on three models: (i) thermonuclear runaway on an accreting white dwarf 
(nova- like event); (ii) He-shell flash in a post asymptotic giant branch star (born-again AGB); and (iii) merger of 
stars. We show that models (i) and (ii) cannot account for the majority of the observed properties of the objects. 
Most significantly, in both nuclear burning type models the object is expected to heat up before declining and 
fade as a very hot compact star. In the observed eruptions the objects declined as very cool giants or supergiants. 
We show that the stellar merger model can account for all the observed properties and conclude that presently 
this is the most promising model to explain the eruptions of the V838 Mon type. 

Key words, stars: supergiants - stars: binary - stars: main sequence - stars: pre-main sequence - stars: mass loss 
- stars: circumstellar matter - stars: variables: general - stars: individual: V838 Mon, V4332 Sgr, M31 RV 



1. Introduction 

The outburst of V838 Mon, discovered at the beginning 
, of January 2002 (Brown 2002}, was an important event 
in the field of stellar eruptions. Initially thought to be 
a nova, the object soon appeared unusual and enigmatic 
in its behaviour. The eruption, as observed in the optical, 
lasted about three months. After developing an A-F super- 
giant spectrum at the optical maximum at the beginning 
of February 2002, the object showed a general tendency 
to evolve to lower effective temperatures. In April 2002 
it almost disappeared from the optical but remained very 
bright in infrared, becoming one of the coolest M-type 
supergiants yet observed. 

The interest raised by V838 Mon recalled two sim- 
ilarly behaving objects observed in the past. These 
were M31 RV, a red variable discovered in M31 in 
September 1988 (Rich et al. 11555) 1 . and V4332 Sgr, dis- 
covered at the end of February 1994 (Hayashi et al. 1994). 
Forgotten for almost 10 years V4332 Sgr appeared, with 
new observations done in 2003, to show very unique spec- 
tral features: a strong emission-line spectrum of very low 
excitation and a huge infrared excess. 



We call M31 RV, V4332 Sgr and V838 Mon V838 Mon 
type objects. Their nature is unclear. As mentioned in a 
number of papers, thermonuclear models (classical nova, 
He-shell flash) seem to be unable to explain these erup- 
tions. This view is, however, not a consensus opinion and 
thermonuclear models are still discussed in the litera- 
ture. Munari et al. J2005) propose a thermonuclear shell 
event in an evolved massive star approaching carbon ig- 
nition, while Lawlor (2005) considers a very late He-shell 
flash in a post- AGB (post asymptotic gaint branch) star 
with a period of accretion added during the flash. Other 
proposed explanations invoke stellar mergers (Soker & 
Tylenda 200.'!; and giants swallowing planets (Retter & 
Marom 125531. 
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In the present paper we make a comprehensive com- 
parison between various models and scenarios proposed 
to explain the nature of the V838 Mon type eruptions 
and the observations of these objects. For this purpose, in 
Sect. |21 we summarize principle observational properties 
of the objects. In Sect. [31 we discuss thermonuclear erup- 
tion models and we show that they cannot consistently 
explain the observations. In Sect. 0] we develop the stellar 
merger model proposed in Soker & Tylenda (2003) and 
show that it can satisfactorily account for all the observed 
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properties of the V838 Mon type objects. A summary is 
given in Sect. 

2. Observed properties 

The principal observational property that makes M31 RV, 
V4332 Sgr and V838 Mon form a separate class of erup- 
tive objects is their spectral and photometric evolution 
in late phases of eruption, in particular during decline. 
The evolution during outburst has been observed only in 
V838 Mon, so we do not know if its light curve prop- 
erties, i.e. multi-outburst character (Crause et al. [2003 
Rctter & Marom 2003,) , are also characteristic of the other 
two objects. Both M31 RV and V4332 Sgr very likely had 
been discovered at a late phase of eruption. Their observed 
evolution was very similar to that of V838 Mon starting 
in last days of March 2002, i.e. about 3 months after its 
discovery. In this phase the objects resambled K-type or 
early M-type supergiants and on a time scale of days they 
evolved to later spectral types (Mould et al. 1990, Martini 
et al. 119991 Tvlenda l2005[l . Soon they became the coolest 
known supergiants or giants. The latest recorded spec- 
trum of M31 RV was M6 1, corresponding to T cff ~ 2600 K 
(Mould et al. H990[) . but after, judging from the subsequent 
infrared photometry, the object evolved to even lower ef- 
fective temperatures. V4332 Sgr reached a spectral type 
as late as M8-9III and a corresponding effective temper- 
ature of ~ 2300 K (Martini et al. [1355) 1 . although it is 
not excluded that the object dropped to even lower T e g. 
In the coolest phase of V838 Mon its spectral type was 
usually being refered to as later than M10 (Desidera & 
Munari 2002) while its effective temperature was as low 
as ~ 2000 K CTvlenda l2UU5}l . 

Before reaching early M-type spectral classes, i.e. for 
T e g <; 3600 K, the objects evolved at a roughly constant 
luminosity. Near the spectral type M0 all the three ob- 
jects began to decline. Initially the luminosity decrease 
was rather fast and a factor of 10 in the luminosity drop 
was reached after ~ 50, ~ 10 and 50-100 days in the 
case of M31 RV (Mould et al. lEOT . V4332 Sgr (Tylenda 
et al. I2005aj> and V838 Mon (TylendaEEJJll, respectively. 
At this luminosity level V838 Mon and presumably also 
V4332 Sgr were close to their minimum effective tempera- 
tures. All that is known about the subsequent evolution of 
M31 RV is that it became fainter by a factor of 100 after 
~ 300 days (Mould et al. HM)| . V838 Mon, after hav- 
ing reached its coolest state, kept declining in luminosity 
but at a much slower rate and with a slowly increasing 
effective temperature. About 1000 days after eruption the 
object was ~ 40 times fainter than at maximum, with 
T e ff ~ 2600 K corresponding to a spectral type of ~ M6 
(TylendalMEl)- At ~ 3500 days after eruption V4332 Sgr 
decreased by a factor of ~ 1500 in luminosity, its spectral 
type was then ~ M3 and T cS ~ 3400 K (Tylenda et al. 
20053- F° r the three objects there were no indications of 
any excursion to high effective temperatures during de- 
cline and no nebular stage was observed. The objects de- 
clined at very low effective temperatures. 



The maximum luminosity attained by an eruptive star 
is an important observational parameter in understanding 
the nature of the object. Unfortunately it depends on the 
distance which is usually uncertain for Galactic objects. 
Therefore M31 RV is the only object for which we can 
be confident of its maximum luminosity of ~ 8 x 10 5 Lq 
(Mould et al. HMty . For a distance to V838 Mon of 8 kpc 
the maximum luminosity of this object was ~ 1.2 x 10 6 Lq 
(Tylenda 12005|. However, because of distance uncertain- 
ties this value is probably uncertain by a factor of ~ 
2. Even more uncertain is the maximum luminosity of 
V4332 Sgr. If its progenitor had been a main sequence star 
the object would have attained ~5x 10 3 Lq (Tylenda et 
al. 12005a)) . which would be significantly lower than in the 
two other cases. 

The outburst amplitude is independent of distance so, 
in principle, it can be determined more precisely than the 
luminosity. Unfortunately only an upper limit is available 
for the brightness of M31 RV before outburst so it can only 
be said that the object brightened by at least ~ 5 mag 
in the / band (Mould et al. 1990). More data are avail- 
able for the progenitors of V4332 Sgr and V838 Mon so 
it can be estimated that at maximum these objects were, 
respectively, ~ 5 x 10 3 and (0.5 — 1.5) x 10 3 times more 
luminous than before the eruption (Tylenda et al. 2005a 
2005b] Tylenda MM- 

The three objects probably lost significant mass. 
However, direct evidence, i.e. P-Cygni line profiles, have 
been observed only for V838 Mon. They give wind ve- 
locities up to ~ 600 km s^ 1 although most of mass loss 
occured at 150 - 400 km s _1 (Munari et al. 12002b! Kolev 
et al. 12*0X721 Crause et al.EUIJJ Kipper et al. "2T7ol)l . Similar 
values, i.e. 270 — 400 km s~ l , result from the observed ex- 
pansion of the effective photosphere in eruption (Tylenda 
2005 ). Widths of emission lines observed in V4332 Sgr, if 
interpreted as primarily due to expansion, indicate a ve- 
locity of ~ 100 km s _1 (Martini et al. I1999|) . Expansion 
velocities of a few hundreds km s^ 1 have been estimated 
for M31 RV (Mould et al lT55f)l Iben & Tutukov HM2l . 

The total mass lost in eruption is uncertain for the 
objects. Very rough estimates done for M31 RV give val- 
ues between 10" 3 - lO^ 1 Mq (Mould et al. [1990). No 
estimate has been done for V4332 Sgr. Rushton et al. 
(2003) obtained an upper limit of ~ 0.01 Mq for the total 
mass lost from V838 Mon from their CO and SiO obser- 
vations. This upper limit should, however, be multiplied 
by at least a factor of 10. Rushton et al. have integrated 
their upper limits for the line fluxes over a velocity in- 
terval of 50 km s , which is much too low compared to 
the observed expansion velocities (see above). Also, their 
estimates assume that all the ejecta were seen in the lines 
(optically thin limit). Tylenda (2005) concludes that the 
total mass lost by V838 Mon is between ~5x 10~ 3 Mq 
and ~ 0.6 M Q . Lynch et al. PJOI) give ~ 0.04 Mq with 
a large uncertainty. An integration of the luminosity de- 
rived by Tylenda (2005) over the pre-eruption and erup- 
tion phase gives a value of ~ 2.5 x 10 46 ergs. Assuming 
equipartition of energy between radiation and mass loss, 
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as well as a wind velocity of 300 km s 1 one obtains mass 
loss of - 0.03 M . 

In the case of V4332 Sgr and V838 Mon inverse P- 
Cygni profiles were observed during the decline, indicat- 
ing an infall velocity of ~ 20 km s _1 (Martini et al. (1999 
Rushton et al. l2005a|l . 

The three objects in decline show strong bands from 
molecules involving oxygen, clearly indicating oxygen-rich 
(C/O < 1) matter. More detailed determinations of the 
chemical composition are available only for V838 Mon 
(Kipper et al. I2UU1 Kaminsky & Pavlenko |2005 ) . The 
abundances of the iron group metals are slightly below 
the solar values, [Fe/H] ~ —0.3. Those of Li, Ba and 
La are somewhat above, although as noted in Kipper 
et al. (|2004ll , the obtained abundances of these elements 
are sensitive to the (uncertain) microturbulcnce velocity. 
Thus it can be concluded that the observed abundances 
in V838 Mon are consistent with unprocessed matter at 
the galactocentric distance of the object. 

Nothing is known about the progenitor of M31 RV. 
In the case of V4332 Sgr, archival photometric data sug- 
gest an F-K type star (Tylenda et al. l2005al Kimeswenger 
2005). V838 Mon has a B3V spectroscopic companion 
(Munari et al. 2002a). The companion is very likely to 
be related to V838 Mon itself, e.g. forming a binary sys- 
tem, although a pure coincidence in the sky cannot be ex- 
cluded. If it is a binary, V838 Mon before the outburst was 
a star similar to its B-type companion, i.e. an early B-type 
main-sequence star or an A-type pre-main-sequence star 
(Tylenda et al. 200§3). The object is surrounded by ex- 
tended circumstellar matter seen in the light echo (Munari 
et al. I2002bl Bond et al. lM^jl . The matter is most likely 
of interstellar origin, possibly being part of one of the star 
forming regions seen in the direction of and at similar dis- 
tances as V838 Mon (TylendaESSl Tylenda et al. l2005b|l . 

Circumstellar matter is also seen in V4332 Sgr, 
strongly radiating in the infrared as observed ~ 9 years af- 
ter outburst (Banerjee et al. lMJSll^UMj) . Most likely this 
matter forms a circumstellar disc dissipating its energy 
in viscous processes (Tylenda et al. 20055} ■ The origin of 
the disc is however unclear. The object also shows an in- 
tense emission-line spectrum due to neutral elements and 
molecules (Banerjee & Ashok l2TJD4l Tylenda et al. l2005a|l . 
A possible interpretation is that this emission is due to 
optically thin (in the continuum) matter left after out- 
burst and now orbiting the central star at a few stellar 
radii (Tylenda et al. I2()05ajl . 

V838 Mon was observed with Chandra a year after 
outburst (Orio et al.|2003 ). The object was not detected 
and an upper limit to the X-ray luminosity was — 0.13 Lq 
for a distance of 8 kpc. 

SiO maser emission was detected in V838 Mon 
— 3 years after its outburst (Deguchi et al. 1 2005 ) . Its in- 
tensity increases with time (Claussen et al. [2005 ). 



3. Problems with thermonuclear eruption models 

A sudden brightening of a star by a factor of — 10 3 is 
usually interpreted as due to a thermonuclear eruption. 
This is because most astrophysical events of this range 
are nova-type outbursts. Their observational properties 
can be consistently explained by thermonuclear runaway 
on an accreting white dwarf. Another mechanism based 
on nuclear burning is a He-shell flash in an evolved star. 
Here we discuss the observed properties of the V838 Mon 
type eruptions in terms of these two thermonuclear mech- 
anisms. 

3.1. Nova-like outbursts 

The V838 Mon type eruptions may be considered as nova- 
like events. The nova mechanism, thermonuclear runaway 
on an accreting white dwarf, has been well studied and is 
known to be able to account for a large variety of astro- 
physical outbursts. This includes relatively "gentle" erup- 
tions lasting decades in some symbiotic stars, slow novae 
that require months to reach maximum and that then can 
stay active for years, and fast novae which in an hour 
or so increase their brightness by 3-4 orders of magni- 
tude and violently eject matter with velocities well above 
1000 km s _1 . Indeed several objects of the nova type, 
like some slow novae, have light curves similar to that 
of V838 Mon. An example is HR Del (Nova Del 1967), 
which showed an initial increase in luminosity that lasted 

— 5 months, then two luminosity peaks with a five month 
separation between them, and then a decline over many 
years (e.g., Drechsel et al. 119771 Rafanelli & Rosino lrF75|l . 
Scaling the luminosity upward and shrinking the time 
scale by a factor of ~ 5 will lead to a light curve qual- 
itatively similar to that of V838 Mon. Friedjung {1992 ) 
pointed out the unusual behavior of this nova in having an 
almost stationary photosphere. This is qualitatively simi- 
lar to V838 Mon, which also seemed to have a quasi-static 
photosphere in its pre-outburst (January 2002) phase 
(Tylenda 2005). However, the light curve is probably the 
only observational feature of the V838 Mon type objects 
that can be easily interpreted in terms of a nova- like event. 
As already noted in early analyses of M31 RV (Mould et al. 
HMIjl and V4332 Sgr (Martini et al.EnHJ), several observa- 
tional properties of these objects are difficult to reconcile 
with the nova hypothesis. 

The crucial argument against the nova-like event 
comes from the observed spectral evolution of the ob- 
jects during outburst and decline. In the case of a classical 
nova event, after an outburst on a degenerate white dwarf, 
steady nuclear burning commences, which halts the lumi- 
nosity typically at — 10 Lq. It can be extinguished only 
after the envelope mass becomes too small to keep the 
burning shell dense and hot enough. But this implies that 
the photospheric radius must be small, hence the effective 
temperature high. This is the reason why novae after erup- 
tion (but also post-AGB stars and late helium shell flash 

- born-again AGB - objects, discussed in Sect. 13. 2|) heat 
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up while evolving at late phases. Models, as well as obser- 
vations, show that T e ff of novae can then reach values well 
above 10 5 K, giving rise to a nebular stage and a signifi- 
cant flux in X-rays. When the nuclear burning stops, the 
object decreases in luminosity as a very hot star. In the H- 
R diagram, the three objects, i.e. RV M31, V4332 Sgr and 
V838 Mon, evolved in the opposite direction to a nova. As 
discussed in Sect. |21 they evolved to lower T e g and faded 
as very cool stars. So far the largest luminosity decrease 
has been recorded for V4332 Sgr, which faded by a fac- 
tor of ~ 1500 in 9 years, still being an M-type object. In 
our opinion, this type of evolution in luminosity and ef- 
fective temperature is evidence that the V838 Mon type 
eruptions are not nove-like events. Below we discuss this 
and other points concerning the nova mechanism in more 
detail and in relation to individual objects. 

Iben & Tutukov (1992 ) have proposed a nova-like sce- 
nario to account for the M31 RV eruption. They consider 
a low-mass (~ 0.6 Mq) cold white dwarf in a short- 
period cataclysmic binary, which accretes at a very low 
rate (~ 10~ n Mq yr _1 ). This allows it to accumulate 
a relatively large mass (~ 0.005 Mq) before thermonu- 
clear runaway. The authors argue that in this case the 
nova outburst would be energetic (maximum luminosity 
well above the Eddington limit) and the object can evolve 
to low effective temperatures due to a massive expanding 
envelope remaining optically thick for a few years. This 
scenario does not allow one to avoid the hot phase after 
the outburst. Sooner or later the expanding envelope be- 
comes transparent and the steady burning white dwarf 
would have to be visible. The phase of steady burning is 
very long for low mass white dwarfs. From the results of 
Yaron et al. (2005) it can be inferred that for an 0.65Mq 
white dwarf it lasts at least ~ 100 yrs. The scenario also 
implies that the binary system is old (~ 10 10 yrs) and that 
the main sequence companion is of low mass 0.3 Mq), 
i.e. an M type dwarf. 

In the case of M31 RV the scenario of Iben & Tutukov 
cannot be conclusively confirmed nor definitively rejected. 
We know nothing about its progenitor. We have no data 
on the object for epochs later than ~ 100 days after its 
discovery in eruption, so a hot phase, when the object is 
expected to be quite faint in optical, cannot, in principle, 
be excluded. However, the hot phase is likely to produce 
strong emission lines when the expanding matter becomes 
ionized. It seems rather unlikely that a nebular stage of 
M31 RV escaped detection in spite of numerous surveys 
for planetary nebulae done for M31 in the last decade. 

In the case of V4332 Sgr the scenario of Iben & 
Tutukov cannot be applied. The progenitor was most 
probably an F-K type star. Thus, if identified as a donor 
secondary, it would imply a system in the upper range 
of the secondary mass (~ 1 Mq) of cataclysmic binaries. 
However, the B — R colour from archive data discussed 
in Tylenda et al. ( 2005a|) is uncertain so an M type pro- 
genitor, required by Iben & Tutukov, cannot be defini- 
tively excluded. But then its absolute magnitude would be 
My <; +9.0, the distance to the object would be 300 pc 



and the luminosity during outburst would be < 150 Lq. 
Such a low luminosity is obviously ruled out in any white 
dwarf thermonuclear event. 

Thus, if continuing with the nova- type scenario in the 
case of V4332 Sgr, an ~ IMq secondary has to be ac- 
cepted. This however implies a typical mass transfer rate 
(3> 10~ n Mq yr _1 ) and a typical nova outburst. Then 
there is no way to evolve to M spectral types after the 
outburst maximum. Second, there is no way to avoid the 
nebular stage. The object is now ~ 1500 times fainter than 
at maximum so the white dwarf had to switch off nuclear 
burning and thus it had to have gone through a very hot 
phase in the past giving rise to a nebular stage. Although 
we have almost no data on the object between 1995- 
2002, it seems rather unlikely that its nebular stage would 
have escaped detection in spite of numerous surveys for 
emission-line objects (e.g. planetary nebulae) done in the 
last decade in the Galactic disc and bulge (galactic coordi- 
nates of V4332 Sgr are I = 13?63, b = -9? 40). A luminos- 
ity drop by a factor of 1000 on a time scale of a few years 
(note that the initial drop in lumonosity by a factor of 
100 occured over 3 months, Martini et al. H999|) requires a 
massive white dwarf in the nova model, which in turn is ex- 
pected to result in expansion velocities ii 1000 km s _1 (e.g. 
Yaron et al. fZQQSI). This is at least an order of magni- 
tude higher than the value of ~ 100 km s _1 observed in 
V4332 Sgr. Finally, the observed radius of the V4332 Sgr 
remnant in 2003, i.e. ~ 5 Rq compared to ~ 1 Rq for the 
progenitor (Tylenda et al. 2005a), implies that the whole 
(hypothetical) nova-like binary would still be embedded in 
a common envelope. Then how was the white dwarf able 
to extinguish nuclear burning if it still has an extended 
hydrogen-rich envelope? 

The progenitor of V838 Mon was a rather hot and lu- 
minous star. If the object forms a (presumably very wide) 
binary with the B3 V companion identified by Munari et 
al. <|2D02a), then the progenitor was either a B2-3 main 
sequence star with a mass ^8 — 10 Mq or a pre-main- 
sequence star of ~ 5 Mq (Tylenda et al. I2005b|) . In the 
former case, one can propose evolutionary paths in bina- 
ries in which a white dwarf can be formed from an ini- 
tially more massive companion, especially if mass transfer 
took place between the binary components. However if the 
cooling time (<; 10 7 years) necessary for the white dwarf 
to cease emitting significant ionizing flux (no significant 
ionization is observed in the diffuse matter surrounding 
V838 Mon) , as well as to become able to produce a nova- 
type outburst, is taken into account, very little room re- 
mains for this possibility. The scenario seems to be very 
unlikely given the fact that we do not know of any cata- 
clysmic variable having a B-type main sequence secondary. 
In the former case the situation is even worse. The age of 
a ~ 5 Mq pre-main-sequence star can be estimated as 
~ 5 x 10 5 years. There is no way to produce and cool a 
white dwarf companion in such a short time. 

If V838 Mon and the B3V companion are observed 
at the same position in the sky due to a random coinci- 
dence then, apart from the above possibilities, the progen- 
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itor of V838 Mon could have been either a BA-type giant 
of mass 2.5 Mq or a B-type post-AGB star (Tylenda 
et al. 2005b). In both cases a white dwarf companion can- 
not be excluded. In the former case the giant would be 
in a rapid expansion phase. Thus if mass transfer begins 
it would increase rapidly to rates close to, or even above, 
the Eddington limit for accretion on a white dwarf. The 
white dwarf would thus be bright and the whole system 
would remind a symbiotic star. In these conditions the 
hydrogen burning on the white dwarf surface would be ig- 
nited rather gently leading to an increase in the system 
brightness only by a factor of a few (compared to the gi- 
ant luminosity plus the white dwarf accretion luminosity). 
V838 Mon increased in luminosity by a factor of ~ 1000. 

The case of a B-type post-AGB star implies that a 
few thousand years ago the object was an AGB giant. 
Thus the only possibility to get a nova-type event is that 
the white dwarf companion was accreting from the AGB 
wind. However, in order to reach a nova peak luminosity 
above 10 5 Lq the accretion rate would have to have been 
<> IQ- 10 Mq yr _1 (Yaron et al. I2UU5)) . Taking a typical 
AGB mass loss rate of 10~ 6 Mq yr _1 (note that at the 
end of AGB, i.e. just before the post-AGB phase, mass 
loss is expected to be significantly higher than that) a 
wind velocity of 10 km s _1 and using Eq. (4.38) in Frank 
et al. (|2002p one finds that the above accretion rate onto 
a 1 Mq white dwarf is possible if the binary separation is 
<; 2 x 10 5 Rq. Such a large binary separation is excluded. 
It would be ii 100 times larger than the maximum effective 
radius of V838 Mon during eruption (Tylenda 2005 ). Thus 
the post-AGB component would be visible all the time, 
including the decline phase. Yet the progenitor was not 
seen in the decline phase when V838 Mon became very 
cool and invisible in the optical (e.g. Munari et al.|J005). 

Thermonuclear outbursts can in principle supply a 
huge amount of energy (e.g., Sparks et al. I1978[l with 
the maximum luminosity at peak of ^ 10 5 Lq achieved 
by fast novae. In V838 Mon the maximum luminosity of 
~ 10 6 Lq lasted ~ 70 days. Novae with a decline of 2 
magnitudes in ti ~ 100 days, have for their maximum vi- 
sual luminosity Ly(max) < 10 5 Lq, and more typically 
V(max) ~ 2 x 10 4 Lq (Cohen HMty . 

The outflow velocity observed in V838 Mon, typically 
~ 150 — 400 km s _1 , can be matched with slow no- 
vae, but not with fast novae. Some models calculated 
by Yaron et al. (2005) reach a nova peak luminosity of 
~ 10 6 Lq. However, the ejection speeds in these cases are 
> 2000 km s _1 , much faster than observed in V838 Mon. 

All the nova models in the extensive grid of Yaron et al. 
(2005) eject < 10~ 3 Mq. In the case of the most luminous 
ones, i.e. having a maximum luminosity ~ 10 6 Lq, it is 
even ^ 5 x 10~ 5 Mq. Estimates for V838 Mon give mass 
loss > 5 x 10~ 3 Mq, most likely ~ 0.01-0.1 Mq. Simliar 
mass loss also probably occured in M31 RV. 

V838 Mon is embedded in diffuse matter which re- 
sulted in the spectacular light echo following the outburst. 
The light echo phenomenon is not a common feature of 
novae. This is what is expected. Fast (previous) ejecta 



should sweep up any circumstcllar or interstellar matter 
from the vicinity of the nova. The light echo has been ob- 
served only in the case of the fast nova GK Per (Nova Per 
1901) (RitcheyEUU Couderc lT355|) . The echoin g dust in 
GK Per was either interstellar in origin (Hessman 1982} or 
due to an old planetary nebula (Bode et al. l2004fl . Tylenda 
( 2004) and Tylenda et al. ( 2005b ) argue that the echoing 
matter in V838 Mon is interstellar. This is compatible with 
the idea that V838 Mon is a young object rather than an 
evolved nova-like system. 

The strong infrared excess observed now in V4332 Sgr 
strongly suggests that the object is surrounded by a large 
(inner radius of ~ 30 central star radii) and luminous 
(~ 15 times brighter than the central star) accretion disc 
(Tylenda et al. 2005a). It is unlikely that a structure like 
this could have formed arround a nova-type binary. 

The abundances usually are important in identify- 
ing the nature of an object. Simulations show that no- 
vae with very large peak luminosity are expected to eject 
highly processed matter with Z up to ~ 0.5 (Yaron et al. 
2005}. As discussed in Sect. [21 the observed abundances 
in V838 Mon are compatible with unprocessed matter. 
Although there are no data on the abundances of key el- 
ements for novae, such as C, N, O, Mg and Ne, no strong 
molecular CN lines were observed in V838 Mon; strong CN 
lines are expected from novae at early stages (Pontefract 
& Rawlings lSUMjl . 

As noted in Orio et al. Q2003[) . the lack of measurable 
X-ray flux in V838 Mon is evidence against a nova-like 
runaway in this object. 

3.2. Helium shell flashes 

While in novae thermonuclear runaway occurs on a rel- 
atively cool white dwarf, in He-shell flashes the runaway 
process takes place in a shell at the outskirts of the core of 
a nuclear active star. This takes place in AGB and post- 
AGB stars. Before the flash the energy generation is main- 
tained by H-burning that adds more helium to the He- 
shell, until conditions for ignition are reached. One of the 
characteristics of this process is that because of expansion 
and cooling, the H-burning is extinguished as the He-shell 
ignites. As a result, the surface luminosity typically varies 
only by a factor of 2-3 during the thermal pulse (Wood & 
Zarro HMTl Vassiliadis & Wood lTMl Her wig 123001) • Thus 
there is no way to obtain an increase of ~ 10 3 in lumi- 
nosity, as observed in V838 Mon and V4332 Sgr, from a 
thermally pulsating AGB star. In the case of V838 Mon 
the AGB hypothesis has to be rejected also on observa- 
tional evidence. As shown in Tylenda et al. (2005b) the 
progenitor of this object was significantly hotter than a 
KM- type giant or supergiant. 

As discussed in Tylenda et al. (2005b) it is possible 
(although unlikely) that the V838 Mon progenitor was a 
B-type post-AGB star. An He-shell flash, usually called 
a late He-shell flash, can also occur at this evolutionary 
stage. However, for the same reason as discussed above 
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(an active H-shell) only a moderate (at most a factor of 
10) increase in brightness is possible (Iben 1984, Blocker 
& Schonberner lT097|) . 

The only way to get a significant luminosity increase 
during a He-shell flash is the so-called very late He-shell 
flash or born-again AGB event which occurs on the cooling 
part of the post- AGB evolutionary track when the H-shell 
is already extinguished (Iben 1984, Herwig 2001, Lawlor 
& McDonald I2003f) . However in this case the object is 
expected to be very hot (~ 10 5 K) before the flash and 
surrounded by a glowing planetary nebula. Indeed, ob- 
jects identified as undergoing a (very) late He-shell flash, 
i.e FG Sge, V4334 Sgr, V605 Aql, are surrounded by plan- 
etary nebulae. Nothing of this kind has been reported for 
the V838 Mon type objects. In the case of V838 Mon itself 
we know from its light echo that the object is surrounded 
by diffuse matter but this matter is not ionized (Orio et al. 
2002J Munari et al. I2002b|) . For the same reason (lack of 
ionization of the light echo matter) the idea of Munari 
et al. Ij2005j) . namely that the eruption of V838 Mon was 
due to a thermonuclear shell flash in an evolved massive 
star (initial mass ~ 65 Mq) close to the Wolf-Rayet phase 
(T e ff ~ 5 x 10 4 K), has to be rejected (Tylenda et al. 
12005 bl> . 

Another important point concerning the late He-shell 
flash is that after the initial excursion to high luminosities 
and low effective temperatures, the object is expected to 
repeat its previous evolutionary path toward high effective 
temperatures. The evolution is thus qualitatively the same 
as that of a nova at late stages. As discussed in Sect. 13.11 
this type of evolution is in a sharp disagreement with the 
observed evolution of the V838 Mon type objects. This 
is one of the crucial arguments against interpreting these 
objects as undergoing a (very) late He-shell flash. 

Finally, as noted in Kipper et al. (20Q3J), the abun- 
dances observed in V838 Mon do not support the idea of 
a late He-shell flash. In particular all the objects identified 
as flashing post- AGB stars (FG Sge, V4334 Sgr, V605 Aql, 
as well as central stars of A 30 and A 78) are C-rich while 
all the V838 Mon type objects have C/0< 1. 

3.2.1. The born-again model of Lawlor 

Lawlor ( 2005 ) proposed a born-again AGB model includ- 
ing an episode of accretion to account for the V838 Mon 
eruption. The accretion is supposed to be due to a burst 
of mass loss from a main-sequence binary companion irra- 
diated by the He-shell flashing white dwarf. The addition 
of accretion produces a second outburst, which the au- 
thor identifies with the main eruption of V838 Mon that 
started in the first days of February 2002. 

The arguments raised above against the very late 
He-shell flash also hold in the case of Lawlor's model. 
This model assumes a very hot progenitor and must 
evolve to high effective temperatures before final fading. 
Other important problems appear when comparing as- 



sumptions and results of the model with the observations 
of V838 Mon. 

As the main observational test Lawlor H2005|) calcu- 
lated a visual light curve from his model and obtained a 
result qualitatively similar to the observed one assuming a 
distance of 6.3 kpc. However when calculating the model 
light curve Lawlor did not take into account interstellar 
extinction which in the case of V838 Mon is as large as 
A v ~ 2.7 (Munari et al. 120051 Tylenda 12005) . Thus an 
agreement between the model and observed light curves 
can be obtained but for a distance as small as 1.8 kpc, 
which can be compared to e.g. a lower limit of 5-6 kpc 
obtained from the light echo expansion (Bond et al. 1^003 
Tylenda 2004, Tylenda et al. 2005bj). For more reason- 
able distance estimates of 8-10 kpc the model of Lawlor 
is systematically too faint by a factor of 20-30. 

Another important problem is that the surface param- 
eters of Lawlor's model just before the accretion event are 
T e g ~ 40 000 K and R ~ 1.5 Rq. This phase is supposed 
to correspond to the pre-eruption phase of V838 Mon 
in January 2002. However, the observed parameters of 
V838 Mon during this phase were T cS ~ 5 000 K and R ~ 
350 Rq (Tylenda 2003J). This discrepancy ensures, as we 
discuss below, that Lawlor's idea of an irradiation-induced 
accretion event cannot occur in the case of V838 Mon. 

Lawlor (2005) proposes that the irradiated compan- 
ion is the B3V companion discovered by Munari et al. 
(2002a). However this companion must be at a distance 
of at least ~ 3000 Rq from V838 Mon (effective radius of 
V838 Mon at the epoch when the companion was discov- 
ered - Tylenda ESQ) and most probably is much father 
away. Obviously, no important irradiation effects nor sig- 
nificant mass exchange can be expected at such a large 
distance between the components. Thus the only possibil- 
ity is that the star identified as the progenitor of V838 Mon 
in Tylenda et al. ( 2005b) was irradiated by its hypotheti- 
cal flashing white-dwarf binary companion. This star most 
probably was an early B-type main sequence star prior to 
the V838 Mon eruption. However, since the irradiation- 
induced mass transfer is supposed to occur during the 
pre-eruption of V838 Mon, i.e. in January 2002, the B-type 
star would have to be at a distance of at least ~ 1600 Rq 
from V838 Mon (supposed to be the flashing white dwarf). 
Otherwise V838 Mon, inflated to ~ 350 Rq in the pre- 
eruption phase, would fill up its Roche lobe and mass 
transfer would occur in the direction opposite to that re- 
quired in the model of Lawlor. Obviously the B-type main 
sequence star would be deeply inside its Roche lobe, so 
in order to obtain significant mass loss from it (Lawlor's 
model requires accretion at a rate of ~ 10~ 3 Mq yr _1 ) 
the irradiated flux would have to be significantly greater 
than the surface flux of the B-type star. This would re- 
quire an irradiation source as luminous as <; 2 x 10 8 Lq. In 
January 2002 V838 Mon had an effective temperature of 
~ 5000 K ( Tylenda 120051 see also e.g. Kimeswenger et al. 
2002, Rushton et al. 12005b)) . i.e. it was much cooler than 
an early B-type star. For obvious reasons a cooler object 
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cannot significantly heat up a hotter one even if the cooler 
one is much more luminous. 

As discussed in Sect.|2V838 Mon lost 0.01 - 0.1 M 
of H-rich matter during its eruption. There seems to be 
no way to get such a large mass loss from the model of 
Lawlor. The mass of the H-rich envelope of a white dwarf 
before a very late He-shell is < 10~ 3 Mq (Iben IT3%l)l 
and < 10~ 4 Mq has been accreted in Lawlor's model. 
Intense mass loss from Lawlor's model would quickly ex- 
pose highly processed matter at the surface. Neither in the 
2002 eruption nor during the decline was there any obser- 
vational indication of this. In particular, the C/O ratio 
remained < 1 all the time. 

We can conclude that the model presented in Lawlor 
(2005) cannot be accepted as an explanation of the 
V838 Mon eruption. It is discrepant with the observations 
in too many points. 

4. The stellar merger model 

In the previous section showed that neither the nova-like 
runaway nor the He-shell flash can account for the ob- 
served properties of the V838 Mon type eruptions. The 
only other way to liberate a significant amount of energy 
in stellar astrophysics is release of gravitational energy, i.e. 
gravitational contraction or accretion. Gravitational con- 
traction powers protostars. It is possible, as discussed in 
Tylenda et al. (|2005bt) . that the progenitor of V838 Mon 
was a pre-main-sequence star. However there is no known 
way in which an energy of <^ 3 x 10 46 erg s~ 4 could be re- 
leased on a time scale of a month in a pre-main-sequence 
star. 

Thus we have to consider accretion events. The more 
compact the accreting object, the more effective is the 
accretion in terms of energy liberated per unit mass ac- 
creted. Accretion onto a neutron star or a black hole, being 
the most effective, is however excluded as it would result 
in an X-ray source. Accretion onto a white dwarf would 
have to precede at a super-Eddington rate to approach 
the observed luminosities and would lead to a nova-type 
thermonuclear event with all the problems involved, as 
discussed in Sect. 13.11 

We are left with accretion on a main sequence star. In 
this case, to obtain a luminosity of ~ 10 6 Lq the accre- 
tion rate must be ^5 0.01 Mq yr _1 . Taking into account 
the total energy involved in the V838 Mon event, as es- 
timated in Sect. 14.11 the main sequence star would have 
to accrete ^0.1 Mq on a time scale of months. An event 
like this is excluded in secular binary system evolution, 
as there is no mechanism that could result in the trans- 
fer of such a large amount of mass in such a short time. 
Young main sequence or pre-main-sequence stars are often 
surrounded by massive protostellar discs. Thermal insta- 
bilities in the disc, similar to those in dwarf novae, can 
transfer a significant mass from the disc to the central 
star. This is the mechanism believed to produce FU Ori 
type eruptions (Hartmann & Kenyon 1996). The discs in 



these objects can transfer as much as 0.01 Mq. However 
this occurs on a time scale of a hundred years. 

Therefore it seems that the only way to get the accre- 
tion event required to explain the V838 Mon type erup- 
tions is accretion of a low mass object onto a main se- 
quence star. This is the reason why we proposed a stellar 
merger model in Soker & Tylenda ( 2003) to explain the 
eruption of V838 Mon. At that time the distance to the 
object was believed to be ~ 1 kpc. Therefore, the merging 
of two low-mass main-sequence stars was considered. 

In the present paper the qualitative general scenario 
for the eruption of V838 Mon remains similar to that in 
Soker & Tylenda (20Q3J. We however change some details 
and update numerical values to be consistent with the 
new distance estimate and more recent data on the evo- 
lution of the object, in particular with the indication that 
the system is more massive and young. A general outline 
of the scenerio presented below (Sect. 14.11) is followed by 
simple simulations of the merger remnant in Sect. 14.21 A 
discussion of the observations within the merger scenario 
is given in Sect. 14.31 

4.1. General considerations and estimates 

We assume, following Tylenda et al. (j2005bf) . that the 
progenitor of V838 Mon was a main sequence star of 
Mi ~ 8A/q or a somewhat less massive star in the pre- 
main-sequence phase. In the following we will call it the 
primary. Most likely it forms a wide binary with a similar 
star, seen as a B3V companion in the decline phase of 
V838 Mon. It is possible that the system also had an- 
other component, or other components, of much lower 
mass. Interactions within the multiple system might have 
dcstabilizated the orbit of one of the low mass compo- 
nents causing it to enter a highly eccentric orbit and, fi- 
nally, to interact with the primary. As the whole system 
is rather young, as discussed in Tylenda et al. ((2005b), 
another possibility is that the V838 Mon progenitor and 
its presumably low mass companion are surrounded by a 
protostellar disc. Interactions between the disc and the 
low mass star might have caused migration of the latter 
and growth of the orbit eccentricity, as discussed, for in- 
stance, in Artymowicz et al. (1991) for stellar binaries or 
in Armitage & Bonnell ( 2002 ) for migrating brown dwarfs 
arround low mass stars, leading to a merger with the pri- 
mary. In the case of a primary at the end of the main 
sequence phase, the onset of fast expansion of the star 
might have contributed to triggering the merger process. 

The mass of the accreted companion, M2, in V838 Mon 
can be estimated from the total energy budget of the 
event. As discussed in Sect. during the main out- 
burst (January - mid- April 2002) V838 Mon lost ~ 2.5 x 
10 46 ergs in radiation and presumably 0.01 — 0.1 Mq in 
mass loss. The energy necessary to lift this matter from 
the surface of a star having Mi = 8 Mq and R\ = 5 Rq 
and to accelerate it to 300 km s _1 at large distances is 
(0.4 - 4.0) x 10 47 ergs. Tylenda pYtB) estimated that 
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the mass of the inflated envelope contracting in the de- 
cline phase was M cnv ~ 0.2 Mq. This value was obtained 
assuming a n = 3/2 polytropic model of the envelope. 
If a n = 3 polytropic model is assumed, as discussed 
below, then the result would be M env ~ 0.1 Mq. The 
energy stored in such an envelope, i.e. a difference be- 
tween the values obtained from Eq. (A. 11) in Tylenda 
lj2"U0l)|> putting i? env = 2000 Rq and R cm = R u is 
(2.5-6.0) x 10 47 ergs. Thus the total energy involved in the 
V838 Mon event is (3 — 10) x 10 47 ergs. Equating this value 
to (GM 1 M 2 )/(2R 1 ) one obtaines M 2 = 0.10 - 0.33 Mq. 

The above analysis can be repeated assuming an 
A-type pre-main-sequence star progenitor discussed in 
Tylenda et al. I|2005b|l . i.e. taking M x = 5 Mq and 
R\ = 7.5 Rq. In this case M env becomes ~ 0.3 Mq and 
~ 0.15 Mq for the n = 1.5 and n = 3 polytropes, re- 
spectively. The mass of the accreted component is then 
obtained as M 2 = 0.15 - 0.5 Mq. 

As discussed in Tylenda et al. (2005a) the case of 
V4332 Sgr can be accounted for by a collision of a planet- 
like component with a ~ 1 Mq main sequence star. 

Let us consider a grazing collision, rather than a 
headon collision, as more plausible in the case of a binary 
merger in an initial eccentric orbit. The gravitational plus 
kinetic energy of the binary system before merger is 



GM1M2 , 

E olh = (1 - e , 

2r n 



(1) 



where r p is the orbital separation at periastron and e is 
the eccentricity. For a grazing merger we have R\ 55 r p ^ 
R\ + R 2 where R2 is the radius of the secondary. In the 
following we will use r p ~ Ri . 

We mostly consider low-mass main-sequence or pre- 
main-sequence stars as secondaries. Their structure can be 
well approximated by a polytropic star with an index n = 
3/2. Thus the total internal (gravitational plus thermal) 
energy of the secondary can be written as (see Eq. IA.2|) 



E-> 



3GMI 
7 R 2 



(2) 



It is reasonable to assume that in the case of a grazing 
collision the low mass secondary, if it is not significantly 
denser than the more massive primary, is disrupted in 
the merger process and forms an inflated envelope of the 
merger remnant. The primary, remaining almost undis- 
turbed, forms a core of the remnant. This scenario is con- 
firmed by the simulations presented in Sect. 14.21 The struc- 
ture of the remnant envelope is determined by the distri- 
bution of the energy dissipated in the merger process. The 
latter is a complex process as demonstrated in numerical 
simulations (e.g. Lombardi et al. 199(3 12002(1 . However, 
to a first approximation we can assume that the rate of 
the energy dissipation is proportional to the density. In 
this case the structure of the remnant envelope can be ap- 
proximated by a uniform energy source model, which is 
close to an n — 3 polytrope (see e.g. Cox & Guili 1984). 
Indeed, the density distribution in the remnant envelopes 
obtained in our simulations discribed in Sect. l4.2l fnot only 



in the cases shown in Fig. ^) can be well reproduced by 
the n = 3 polytropic distributions (see Eq. A4 in Tylenda 
2005). For an outer radius of the envelope, i? C nv, much 
larger than the inner one (assumed to be the radius of the 
primary, i?i) the total energy of the envelope, E env , can 
be approximated by (see Eq. A. 17 in Tylenda 
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2R 1 ln(R em /R 1 ) 
The energy balance can be written as 
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which, using Eqs. ~ ©, can be transformed to 
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In the above energy balance we do not consider energy 
losses in form of radiation and mass loss. In the global 
energetics of the event they can be important, as shown in 
our estimates above for V838 Mon. However, these losses 
mostly take place in later phases of the evolution of the 
merger remnant, which occur on a thermal time scale. The 
merger is expected to occur on a dynamical time scale so 
for the initial structure of the merger remnant the losses 
are not expected to be important. This is consistent with 
numerical simulations of stellar mergers (e.g. Lombardi 
et al.Enni see also Sect. l4~2"j) . 

Analyses of observations give R cnv /Ri ^> 1. Maximum 
values of the photospheric radius in the case of V838 Mon 
reached ~ 3000 Rq (Tylenda EESt This however refers 
to freely expanding matter rather than to a hydrostatic 
remnant. The latter was smaller but probably still as large 
as ~ 2000 Rq, as inferred from the photospheric radius 
observed during the early decline. For R\ = 5 Rq this 
corresponds to R cnv /Ri ~ 400. In the case of V4332 Sgr 
the envelope was probably less inflated but still as large 
as R cnv /Ri ~ 100 - 150 (Tylenda et al. HW05ajl . 

Eq. J5J) shows that from a merger of two similar stars, 
e.g. two main sequence stars of not vary different mass, 
we do not expect a significantly inflated remnant. For 
instance, taking Mi — M 2 and R\ ~ R 2l Eq. JSJ) pre- 
dicts 1.7 ^ R cnv /Ri 3.2 for < e < 1. Note, however, 
that the cases of two similar stars are not particularly ap- 
plicable to the above approach, which assumes an undis- 
turbed primary and a largely inflated remnant envelope. 
Simulations, described in Sect.Olfor two 0.6 Mq ZAMS 
stars give 6.5 < R cnv /Ri < 8.4 for 0.5 < r p /(Ri + R 2 ) < 
1.0, thus the remnants are more extended than predicted 
from Eq. JSJ but still much less than observed. 

Taking M x = 8Mq, M 2 = 0.3(0.1) Mq, i.e. masses 
expected for V838 Mon, and the corresponding main se- 
quence radii, i.e. R\ = 5 Rq and R 2 — 0.35(0.15) Rq, 
one still gets little expansion from Eq. JSJ, i.e. R cnv /Ri ^ 
9.0(16.5). However the above set of stellar parameters is 
not realistic. From evolutionary considerations it can be 
concluded that at an age of an 8 Mq star leaving the main 
sequence, stars of masses ^1.0 Mq still have to be in the 
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pre- main-sequence phase. As can be seen from Eq. Ip|l. 
larger values of R 2 would result in more inflated merger 
remnants. An 8Mq star reachs the ZAMS at an age of 
3.1 x 10 5 years (see Sect. IP)) . At this age a 0.3(0.1) M 
pre-main-sequence star is expected to have a radius of 
~ 2.5(1.3) Rq, as estimated from Eq. (|A.4|I . To get enve- 
lope extensions as large as R env /Ri ~ 400 Eq. J5j requires 
i? 2 ^ 0.95(0.32) Rq. 

If the V838 Mon progenitor was an early A-type pre- 
main-sequence star, we can take Mi ~ 5 Mq and R\ ~ 
7.5 Rq. The age of the star would be ~ 5 x 10 5 years 
(see Sect.E21>. At this age the radius of a 0.5(0.15) Mq 
pre-main-sequence star would be ~ 3.5(1.5) Rq. For 
these parameters (and assuming e = 1) Eq. JSJ predicts 
i? onv /i?i ~ 230(2400). 

In the case of V4332 Sgr we can assume M\ = 
IMq, Ri — 1 Rq and R2 = 0.1 Rq (radius typical 
for brown dwarfs and Jupiter- like planets). Then infla- 
tions of R cnv /Ri ~ 150 can be obtained from Eq. (j^J if 
M 2 < 0.02 M . 

Even in the most favorable cases, i.e. when the term 
(M 2 /M 1 )(R 1 /R 2 ) in Eq. © is negligible, e ^ 0.8 is re- 
quired to get the remnant extension as large as R cnv /Ri iS 
150. 



4.2. Numerical simulations 

We have carried out simulations of a merger of 
two stars with parameters expected for V838 Mon. 
We used the MMAS (version 1.6) package de- 
scribed in Lombari et al. I|2002l I2()()3|l and available 
at http:/faculty.vassar.edu/lombardi/mmas/. This code, 
based on simple physical arguments and algorithms, al- 
lows one to obtained a one dimensional structure of the 
collision product without hydrodynamic simulations. It 
still accounts for shock heating, mass loss and fluid mix- 
ing. The remnants returned by MMAS are in hydrostatic 
equilibrium but before thermal relaxation. The code ex- 
plicitely assumes e = 1.0 before collision so its results can 
be considered as reasonable for highly eccentric binary or- 
bits. 

Following the discussion in Sect. 14.11 we have consid- 
ered an 8.0 Mq main sequence primary which accretes 
an 0.3 Mq pre-main-sequence secondary. As argued in 
Sect. 14.11 a grazing merger is more likely than a headon 
collision, so we discuss cases with 0.5 < x p < 1.0, where 
x P = r v /{R\ + R 2 ). 

The structure of the primary before merger 
was obtained from the TYCHO (version 6.0) stel- 
lar evolution code (Young et al. 12001(1 available at 
|http:/ /chandra.as.arizona.edu/^dave/tycho-intr o.htmH 
Starting from an initial (pre-main-sequence) model 
the 8.0 Mq star reaches the ZAMS at an age of 
3.15 x 10 5 years. For our simulations we have taken the 
stellar structure at 5.7 x 10 6 years. The star is well on 
the main sequence with ~ 10% of hydrogen burnt in the 
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Fig. 1. Results of a grazing collision of two stars using the 
MMAS package. The initial system consists of an 8.0 Mq 
main sequence primary and an 0.3 Mq pre-main sequence 
secondary at an evolutionary age of 5.7x 10 6 years. Dashed 
curves: density profiles of the stars (upper: primary - 
lower: secondary) before collision. Full curves: density pro- 
file in the merger product. Dotted curves: values of a pa- 
rameter showing the origin of the matter in the merger 
product (—8: matter from the primary, —12: matter from 
the secondary - see text). The curves (full and dotted) are 



marked with the value of ; 



core and an outer radius Ri = 3.7 R 



©• 



The structure of the secondary has been approximated 
by a n = 3/2 polytrope, as discussed in ApnendixlAl Using 
Eq. (|A.4J| one gets that at an age of 5.7 x 10 6 years the 
radius of the 0.3 A/q secondary is R 2 — 0.9 Rq. 

Density profiles for the stars (before merger) modelled 
as explained above are shown with dashed curves in Fig.^ 
Full curves present the density profile of the merger prod- 
ucts obtained with MMAS. Three cases of x p = 0.6, 0.8 
and 1.0 are shown. We have defined a parameter to which 
values of —8 and —12 have been assigned for the matter 
in the primary and secondary, respectively, before merger. 
MMAS sorts and mixes the matter from both components 
after the merger event so the returned value of this param- 
eter allows us to see from which component and in what 
proportion the matter comes at a given radius of the rem- 
nant. This result is shown with dotted curves in Fig. ^ 

In the case of a main sequence star the density gradient 
is relatively flat in the interior, that quickly drops when 
approaching the star surface. Therefore if the star suffers a 
grazing collision with a low mass pre-main sequence star 
its interior structure remains almost intact (see Fig. 0. 
Only outer regions are seriously disturbed and together 
with the matter from the destroyed secondary they form 
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an inflated envelope. 1 The outer radius of the envelope 
is very sensitive to x p and it becomes as large as 3.7 x 
10 3 Rq in the case of x p = 1.0 shown in Fig. 2] Note 
that Eq. (J5J predicts R cnv ~ 7 x 10 3 Rq which is in a fair 
agreement with the simulations. The mass of the inflated 
envelope is ~ 0.30 Mq. For x p ^ 0.8 the remnant envelope 
is predominantly formed from the disrupted secondary. 
The simulations predict ~ 0.015 Mq lost from the system, 
for x p <; 0.8 almost entirely from the secondary. 

We have also carried out simulations for the same stars 
as in Fig. ^ but at an age of 3 x 10 7 years as well as for 
two pre-main-sequence stars of 5.0 and 0.3 Mq at an age 
of 5 x 10 5 years. In the former case the 8 Mq primary just 
started to move from the main sequence and the radii of 
the components were R\ — 9 Rq and R^ — 0.6 Rq. In the 
latter case the effective temperature of the 5 Mq primary 
was ~ 10 000 K, i.e. corresponding to an ~ AO spectral 
type discussed in Tylenda et al. (2005b). The stellar radii 
were R± = 8 Rq and R2 = 2 Rq . The results were very 
similar to those displayed in Fig. ^ In particular the radii 
of the obtained remnants were 10 2 — 10 4 Rq . 

In all the cases discussed the density distribution in the 
inflated envelope of the remnant can be well approximated 
by a polytropic envelope (see Appendix in Tylenda[2005) 
with an index n = 3. 



4.3. Discussion 

As we have shown using simple estimates in Sect. I4~T1 and 
simulations in Sect. 14.21 it is possible to obtain merger 
remnants with radii similar to those observed in V838 Mon 
and V4332 Sgr. These results should be verified with more 
sophysticated modelling. As noted in Sect. l4.2l thc MMAS 
code is a simple algorithm. It has been calibrated on SPH 
simulations of mergers of two low-mass (0.4 — 0.8 Mq) 
main-sequence stars. For the star parameters used in 
Sect. 14.21 the code works at its limits and the results are 
probably subject to considerable uncertainties. 

Some important processes are not taken into account 
in the MMAS code. These include effects of radiation 
transfer and angular momentum. Radiation pressure be- 
comes particularly important if the luminosity is super- 
Eddington. This was certainly the case in V838 Mon. In 
February-March 2002 the observed luminosity was above 
the Eddington value (Tylenda 12005(1 . which means that 
the rate of energy dissipation deep in the envelope must 
have been super-Eddington by a larger factor. A part of 
energy went to expansion of the envelope and mass loss. 
Also energy diffusion processes in the envelope work in the 
sense that the observed eruption lasts longer, while the 
maximum luminosity is lower, compared to the original 
energy burst. Indeed, an analysis done in Tylenda (p005) 
suggests that the main eruption observed in February- 



1 The situation is qualitatively different for collisions ap- 
proaching headon, i.e. x p 0.4. Then most of the secondary 
dives into the primary centre and forms the core of the collision 
product. 



March 2002 was generated on a time scale of a few days 
at the end of January 2002. 

The rate of energy dissipation in a merger event can 
be estimated from 
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p V 3 7T i?2, 



(6) 



where v is the velocity of the accreted component of radius 
i?2 moving in an ambient medium of density p. Assuming 
that the merger takes place in the outer layers of the pri- 
mary and that v is close to the Keplerian velocity, i.e. 
v ~ yjGMi/Rx, Eq. © becomes 
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This is a very rough estimate but shows that in the case 
of mergers simulated in Sect. 14.21 the energy is expected 
to be dissipated at a very high rate, orders of magnitude 
above the Eddington limit. This can have important ef- 
fects on the structure of the remnant and mass loss. In 
particular the total mass lost in the merger can be signif- 
icantly higher than the values predicted from dynamical 
effects only in the MMAS code (0.01 - 0.02 Mq in our 
cases). 

The evolution of a merger remnant can be quite easily 
predicted, at least qualitatively. The high energy sources 
will quickly decline after the merger. The rate of the en- 
ergy transported from the base of the inflated envelope 
will soon become limited to the luminosity of the central 
star, which is much less than the luminosity radiated away 
from the photosphere. As a result, the matter in the en- 
velope, which did not achieve escape velocity during the 
merger eruption, has to contract under the gravity of the 
central star. In this way the internal energy of the con- 
tracting envelope can be liberated at a rate necessary to 
balance the photosphcric energy loss. 

A transition between the eruption phase (dominated 
by expansion of the matter) and the decline phase (domi- 
nated by gravitational contraction of the envelope) is ex- 
pected to occur on a dynamical time scale, which is 



R : : 



2 G Mi 



ill 



(8) 



For the observed radius of V838 Mon at the beginning of 
the decline, i? e nv — 2000 Rq, (and assuming Mi = 8 Mq) 
Eq. JSJ) gives td — 1 year. Taking i? cnv ~ 150 Rq and 
Mi = 1 Mq as the parameters for V4332 Sgr (Tylenda 
et al. I2005a[) one gets — 1 month. During this transi- 
tion phase the photospheric regions are deprived of energy 
supply (merger energy is declining, contraction not yet ef- 
fectively started) so they must cool down due to radia- 
tive processes. In V838 Mon this phase was observed over 
~ 200 days since ~ 10 April 2002. As can be seen from 
Fig. 2 in Tylenda |2DD2J, the object kept the effective ra- 
dius between 2000 — 3000 Rq, its effective temperature 
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dropped from ~ 3600 K to a minimum value of ~ 1800 K, 
while the luminosity declined by more than an order of 
magnitude. In the case of V4332 Sgr the transition phase 
probably started on - 5 April 1994 (Tylenda et al. l2005a|) . 
During the subsequent three months the object dropped 
by a factor of ~ 100 in luminosity, reaching an effective 
temperature as low as ~ 2300 K (Martini et al.^999). 

In later phases the evolution of a merger remnant will 
be governed by the gravitational contraction of its inflated 
envelope. Thus the main processes are here the same as in 
protostars. As a result we expect that on the HR diagram 
the object should decline along a Hayashi track. However 
the time scale of evolution of the merger remnant must 
be much shorter than that of a protostar, as instead of 
the whole star, only a small mass envelope is involved in 
the contraction. As shown in Tylenda l|2UU5fl and Tylenda 
et al. J2Q05a) the observed evolution of V838 Mon for ~ 
2.5 years after its 2002 eruption and the observed state of 
V4332 Sgr ~ 9 years after its 1994 eruption are consistent 
with this prediction. 

This behaviour of V838 Mon type objects in initial 
decline resemble the expected evolution of supernova type 
la companions as discussed in Podsiadlowski f2j0n3$. These 
stars, when puffed up to large radii and luminosities by 
the supernova explosion, due to a very short thermal time 
scale are expected to enter into the "forbidden" Hayashi 
region and then decline rapidly. 

Angular momentum is not taken into account when 
determining the remnant structure in the MMAS code. 
The code however provides the expected angular momen- 
tum distribution in the remnant. In our cases the ratio 
of the resultant specific angular momentum to the local 
Keplerian value is roughly constant throughout the in- 
flated envelope and is typically 0.10-0.15. This may have 
important effects on the structure of the remnant, partic- 
ulary on its evolution with time. The above result means 
differential rotation of the envelope with the rotation ve- 
locity increasing inward as ~ r -1 / 2 . Radial motions, for 
instance convection, would transport angular momentum 
outward, thus accelerating rotation in outer layers. This 
may lead to an expansion of the envelope and/or intensi- 
fied winds from the remnant in the equatorial directions. 
In later epochs, when the envelope as a whole is expected 
to contract, the angular momentum may prevent equato- 
rial regions from significant contraction. Thus some matter 
can be left at larger distances, forming a Keplerian ring or 
disc arround the contracting central object. As discussed 
in Tylenda et al. (|2005a|l this is a likely explanation of the 
origin of cold matter seen in the emission molecular bands 
and atomic lines in the spectrum of V4332 Sgr. 

One of the principal observational characteristics of an 
eruptive object is its light curve. V838 Mon displayed a 
quite complex light curve in January- April 2002, which 
can be interpreted as a sequence of three consecutive out- 
bursts. This led Retter & Marom l|2003|l to propose that 
the V838 Mon eruption was due to a red giant that swal- 
lowed three planets. This idea however suffers from sev- 
eral important shortcomings. First, as shown in Tylenda 



et al. (|2005b|l the progenitor of V838 Mon was not a red 
giant. Second, the energy budget of the V838 Mon erup- 
tion in Sect. 14.11 is such that a mass two orders of mag- 
nitude larger than three Jupiter-like planets is required 
to account for it. Third, it seems impossible to propose a 
scenario in which three planets would fall onto the cen- 
tral star one by one on a time scale of months. A study 
in Tylenda (2005) shows that if a difference in propaga- 
tion time between the second (observed in February 2002) 
and the third (observed in March 2002) outburst is taken 
into account then the conclusion is that the two out- 
bursts were generated on a time span of a few days only. 
The whole outburst can be divided into two main phases. 
The first one, observed in January 2002 and called pre- 
emption in Tylenda (2005), was marked by a relatively 
static photosphere of ~ 350 Rq. The second one was ob- 
served in February- April 2002, called eruption, and was 
marked by an expansion of the photosphere from ~ 350 
to - 3000 Rq. 

The question that arises is how this multi-outburst 
light curve of V838 Mon can be explained in our merger 
model. Two general interpretations can be proposed to 
account for it. The first one is that the event involved 
accretion of two companions on a time span of a month. 
The second is that it was a single companion that collided 
and merged with the V838 Mon progenitor but that the 
merger event proceeded in multiple phases. 

Accretion of two companions can occur in two scenar- 
ios. One possibility is that orbits of two low-mass compan- 
ions had been disturbed, say due to a destabilizing event 
in the V838 Mon system, in such a way that the two com- 
panions collided with the primary star one after the other. 
Another possibility is that only one companion had been 
disturbed in its orbit so it collided with the primary but 
the merger event dcstibilizcd another companion causing 
its merger with the primary a month after the first one. 

A more detailed scenario in the former case is as fol- 
lows. The first companion collided with the primary star 
at the end of December 2001 and initiated the pre-eruption 
phase. As a results an envelope inflated to ~ 350 Rq was 
formed. If there were another companion orbiting the pri- 
mary at a distance much smaller than the radius of the 
envelope then it would become engulfed by the envelope. 
Drag forces between the orbiting companion and the enve- 
lope caused the companion to spiral in, so after a month it 
reached dense regions of the primary, more or less undis- 
turbed by the accretion of the first companion. Then the 
second companion was disrupted, liberating its orbital en- 
ergy which resulted in the main eruption started at the 
begining of February 2002. However, as we show below, it 
is very difficult, perhaps even impossible, to get a spiraling 
time as short as 1 month for realistic parameters in this 
scenario. 

The luminosity of V838 Mon in the main eruption was 
~ 20 times higher than that in the pre-eruption. This im- 
plies that the second companion was more massive than 
the first one by a factor similar to the luminosity ratio. 
Thus we can assume that the mass of the first compan- 
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ion was ~ 0.01 Mq. Let us approximate the structure of 
the envelope resulting from its disruption by a polytropic 
model (see Appendix in Tylenda 2005 ) . The radius of the 
envelope was ~ 350 Rq, as observed in the pre-eruption 
phase, while its mass can be assumed to be comparable to 
the mass of the distrupted companion. Taking the mass 
and radius of the primary, Mi = 8 Mq and R\ — 5 Rq, 
we can estimate a density at the base of the envelope, 
po, from Eq. (A. 5) in Tylenda ( 2 0T)5jl . In the case of the 
polytropic index n = 3/2 one gets pa ~ 3 x 10~ 7 g cm -3 , 
while for n — 3 the result is po ~ 1.5 x 10~ 5 g cm~ 3 . 
Now we can use Eq. (|B.7|) to estimate the time neces- 
sary for the second companion to spiral in, t s . Assuming 
M 2 ~ 0.3 Mq, R 2 ~ 0.9 Rq (as in Sect.OJl and taking a 
rather lower limit to the orbital separation of ~ 2i£i, i.e. 
putting r /R a ~ 2.0 in Eq. (|R7|) . the result is t s ~ 150 
and 6 years for n = 3/2 and 3, respectively. As can be seen 
from Eq. i|B.7|l t s strongly increases with increasing values 
of ro/Ro, e.g. for ro/Ro = 10 the in-spiral time becomes 
5000 and 2000 years, respectively. 

Another point is that during the spiraling phase the 
companion releases energy ~ [GMiM 2 (l — i?o/?"o)]/ (2i?o) 
which for i?oAo -C 1 becomes comparable to the energy 
released during the disruption of the companion at Rq. 
Thus even if it were possible to spiral in the companion 
over a month, the energy released during the pre-eruption 
phase would be comparable to that in the eruption phase. 
In the case of V838 Mon the ratio of the energy radiated 
away during the pre-eruption phase (January 2002) to that 
in the eruption phase (February - mid-April 2002) was 
- 0.03. 

Therefore we can conclude that if the 2002 eruption of 
V838 Mon was due to accretion of two star-like objects the 
only possibility seems to be that the two bodies arrived 
from outside on very similar trajectories: the lower-mass 
one was followed by the more-massive one with a month 
between them. It could have been a low-mass binary (e.g. 
a ~ 0.3 Mq star orbited by a ~ 0.01 Mq brown dwarf) 
"injected" by dynamical effects from the periphery of the 
V838 Mon system or from outside it. 

In our opinion, the second interpretation is more rea- 
sonable. The outburst was likely due to the merger of a 
single companion but that the merger proceeded in mul- 
tiple phases. In several numerical simulations of star colli- 
sions described in the literature this is indeed the case (e.g. 
Lombardi et al. |T996, Freitag & Benz 2005 , see also ani- 
mations available at http://www.manybody.org/modest). 
The components, often partly disrupted after the first en- 
counter, leave each other on a highly eccentric orbit and 
collide again at the return. This may be repeated, some- 
times several times, with a rapidly decreasing time be- 
tween the encounters. Thus a plausible interpretion of the 
light curve of V838 Mon is as follows. The first major en- 
counter betweeen the low-mass companion and the proge- 
niotor of V838 Mon happened at the end of December 2001 
and produced the pre-eruption phase. The companion 
probably became partly disrupted and only a small part 
of the available energy was dissipated. Matter significantly 



shocked during the first encounter formed the pre-eruption 
envelope. Most of the matter from the companion collided 
with the primary after a month (note that the dynam- 
ical time scale - Eq. (jSJ - for the observed size of the 
envelope in the pre-eruption is ~ 30 days) and once again 
after a few days. This produced the main eruption with 
the double-maximum structure in the light curve observed 
in February and March 2002 and two shells in mass loss 
found by Tylenda pTHS)) . 

Before the main encounters, which led to the observed 
outburst, it is possible that there had been a series of 
minor encounters due to the companion passing by the 
periastron. This could have led to minor outbursts and a 
long-term increase in brightness before the main outburst. 
In the case of V838 Mon we have no observational data be- 
tween 1994 and the outburst discovered at the beginning 
of 2002. Before 1994 the object had remained constant, 
at least for ~55 years (Goranskij et al. 2004). In the case 
of V4332 Sgr, however, Kimeswenger f2005|) shows that 
the object started to rise several years before the outburst 
observed in 1994. This can be interpreted as a result of 
minor interactions before the final merger. If V4332 Sgr 
is a Galactic buldge object, then its progenitor was a G- 
type giant, as discussed in Tylenda et al (2005a). Then 
the merger process may have occured well below the gi- 
ant photosphere. Before the merger, when the companion 
was spiraling in the giant envelope, the brightness of the 
object was expected to steadily increase, as can be seen 
from Eq. I|B.9|I . 

We now discuss the observed abundances in the 
V838 Mon type objects. As discussed in Sect. there 
is no observational indication that the matter observed in 
the objects was significantly processed by nuclear burn- 
ing. In our merger model the observed matter is primarily 
due to the disrupted low-mass secondary. In the case of 
V838 Mon we argue that it was a low-mass pre-main- 
sequence star so no significant nuclear processing took 
place inside it. In particular this concerns the observed 
abundance of Li. Low-mass stars of masses ^ 0.5 Mq are 
not expected to have significantly burnt Li at an age of 
^ 10 7 years (Chabrier & Baraffe H997l also our own mod- 
elling done with the TYCHO code). Even if some mass 
from the primary had been mixed into the remnant it 
should have no significant effect on the surface Li abun- 
dances. Li is expected to remain unburnt in outer layers 
of an 8 Mq main-sequence or a 5 Mq pre-main-sequence 
star. 

5. Summary and discussion 

Table ^ summarizes our discussion of the models com- 
pared with the observations. For clarity it is restricted to 
the observed properties of V838 Mon only. Column (1) 
lists the principal observational properties of the object. 
Comments on how these properties can be explained in the 
classical nova, born-again AGB (very late He-shell flash) 
and stellar merger models are given in columns (2) - (4), 
respectively. In the table we do not discuss the idea of 
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Observed property 




Classical nova 


Born-again AGB 


Stellar merger 


(1) 




(2) 


(3) 


(4) 


Increase by factor ~ 


10 3 


(+) a Possible, although 


(+) Possible if starts well 


(+) Easy to obtain 


in luminosity (AV — 


:7.5) 


usually larger 


down the cooling track 




Multi-outburst 




(+) Compatible with 


( — ) Cannot be explained 


(+) Multiple phase 


light curve 




a slow nova 


(Lawlor's accretion episode 
cannot work) 


merger process 


Fading as a very 




( — ) Contrary to expected 


(— ) Contrary to expected 


(+) Contraction along 


cool supergiant 




(decline at T cS Z 10 5 K) 


(decline at T cff Z 10 s K) 


the Hayashi track 


L ~ 10 b L Q 




( — ) Too luminous for a slow 


( — ) Too luminous 


(+) Merger of ~ 0.3M Q 


over ~ 70 days 




nova, too long for a fast nova 


and too fast 


and ~ 8Mq stars 


^ solar abundances 




(— ) ~ 10 x solar abundances 


( — ) Enrichment by 


(+) Young stars at 


([Fe/H] ~ -0.3) 




expected 


dredge-up expected 


galactic outskirts 


Outflow velocity 




(+) Observed in slow novae 


(-) Too fast 


(+) ~ escape velocity 


200 - 500 km s" 1 








from a MS star 


B-type progenitor 




(— ) Cannot work 


(— ) Too cool for a star 
before born-again AGB 


(+) B-type MS primary 


Association with 




( — ) too little time to form 


( — ) Incompatible with 


(+) Progenitor of similar type 


a young (B3 V) star 




and cool a WD 


a post- AGB progenitor 




Circumstellar 




( — ) Past eruptions expected 


(-) Past AGB wind 


(+) Part of an ISM region 


non-ionized matter 




to sweep up the enviroment 


expected to be ionized 





a (+) can be explained by the model; (— ) cannot (very diffi- 
cult to) be explained by the model 



an evolved very massive star proposed in Munari et al. 
(2005). This is mainly because this idea is ruled out by 
the observed lack of noticeable ionization of the matter 
seen m the light echo (Tylenda et al. l2005b|) . There is no 
available model based on this idea so it is not clear what 
observational properties are predicted from it. 

Tabled shows that presently the merger model is the 
most promising one for V838 Mon and the other objects 
of the same type. Both classical nova and He-shell flash 
models seem to be ruled out by too many contradictions 
with the observed properties of the objects. 

The main drawback of our merger model is that it is 
based on estimates and approximate considerations. The 
simple simulations presented in Sect. 14.21 only demon- 
strate that an inflated envelope can be formed from 
a merger. Realistic simulations should provide a three- 
dimensional structure of the inflated envelope, its evo- 
lution after merger, and the expected light curve. This 
would however be a very complex task involving three- 
dimensional radiation- hydrodynamics. We hope that this 
paper will stimulate research in this direction. 

So far we have seen three events of the V838 Mon 
type: two in our Galaxy, one in M31. A binary population 
synthesis done in Hen, Podsiadlowski & Eggleton (1995) 
shows that 1-2 stellar mergers are expected in our Galaxy 
per 10 years. This result considers only binaries with stel- 
lar components having masses > 0.8 Mq. Thus if multi- 
ple systems, lower mass stars, brown dwarfs and massive 
planets were also considered, two events in our Galaxy in 
a time span of 8 years are not surprising. What is supris- 
ing is that no event of this kind had been reported (in our 
Galaxy) before the discovery of V4332 Sgr. Kato (2003) 



suggests that Nova CK Vul 1670 might have been a stellar 
merger event. Bond & Siegel (|2005[) mention that Nova 
V1148 Sgr 1943 had a late-type spectrum thus being a 
possible object of the V838 Mon type. It is quite probable 
that some V838 Mon type objects remain hidden among 
stars catalogued as classical nova or nova-like. 

Mergers of massive stars are considered as a possi- 
ble way to form very high mass stars. Bally & Zinnecker 
( 2005]) have discussed expected observational appearances 
of mergers involving massive stars. These events can re- 
lease up to 10 51 ergs and reach a peak luminosity above 
10 7 L Q . They argue that an event of this kind happened 
~ 500 years ago in the Orion molecular cloud. Their study 
as well as our work shows that stellar mergers are impor- 
tant not only to understand the evolution of binaries or 
the origin of blue stragglers. Due to the large amount of 
energy released on a short time scale they can also pro- 
duce spectacular events, easy to observe even beyond our 
Galaxy. 

Appendix A: Low mass pre-main-sequence stars 

Low mass pre-main-sequence stars are convective. In this 
case, for given values of the star mass, M», and radius, 
the stellar structure can be obtained from a solution of the 
Lane-Emden equation with a polytropic index n = 3/2 
(for detailes see e.g. Cox & Giuli 11984")) . 

A pre-main-sequence star radiates mainly due to its 
gravitational contraction, so its lifetime, i pms , can be es- 
timated from (see e.g. Stahler 1994) 

— £* 

^pms — — j • (A-1) 
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Here E* is the total (gravitational plus internal) energy of 
the star, which, for an n = 3/2 polytrope, is (see e.g. Cox 
& GiuliHHEH> 



3 GMi 



7 R* 1 
while the luminosity, L*, is 

L. =\vRloTt, 



(A.2) 
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where T* is the stellar effective temperature. Combining 
the above equations one gets 
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where / pms accounts for an uncertainty due to Eq. jA.ljl . 
Comparing results of evolutionary model calculations for 
low mass pre-main-sequence stars (StahlerHnni our own 
models calculated with TYCHO) with the predictions 
from Eq. (|A.4() we have found / pms = 0.6 - 0.8. 

Appendix B: Spiraling in a polytropic envelope 

Let us consider a companion of mass M 2 and radius R2 
that orbits inside a non-rotating envelope of a primary 
star of mass, M\, at a distance, r, from the star centre. 
We assume that M 2 <C M\ and that the density in the 
envelope, p, is much lower than the mean density of the 
companion, p 2 = {3M 2 )/{&-kR\). Thus we can neglect ef- 
fects of accretion or evapration processes on the mass of 
the companion. Then we can assume that the compan- 
ion moves with a Keplerian velocity, v — yj GMi/r, so its 
energy (gravitational plus kinetic) is 

G Mi M 2 



E 



2r 



(B.l) 



Due to drag forces the companion loses its energy at a rate 
which we approximate by 



— — ~ — ir R 7 pv — . 
dt 2F 2 
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Taking a time derivative of Eq. (|B.1|I and substituting it 
in Eq. I)B.2|1 one gets the spiraling rate 
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We assume that the structure of the envelope is given 
by a polytropic solution of index n so, apart from the 
outermost regions of the envelope, the density distribution 
can be approximated by (see Eq. A.4 in Tylenda 
for x <C 1 and xq <C 1) 



(Ra 
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where po is the density at the base of the envelope at 
r = Rq. Then Eq. 1B.3() can be rewritten as 



- 1 ' 2 dr 
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which can be integrated giving 

1\ vr R% 



r " +1/2 - r " +1 /2 = [n +- 



M 2 



p R™ {G Mi) 1 / 2 1, (B.6) 



where tq is the radius of the orbit of the spiraling com- 
panion at t = 0. Eq. (|B.6j) can be used to obtain the time, 
t s , necessary for the comapnion to spiral from the initial 
orbit at r$ to the base of the envelope at r = Rq, i.e. 
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where 
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where Po = 2ttRq tJRq/(GM\) is the orbital period at the 
base of the envelope. 

Eqs. (|B.6(I and Eq. l|B.2j) can be used to obtain the 
luminosity generated by the spiriling companion 
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